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ABSTRACT 

We present a sample of 131 quasars from the Sloan Digital Sky Survey at redshifts 0.8 < z < 1.6 
with double peaks in either of the high- ionization narrow emission lines [Ne V]A3426 or [Ne III] A3869. 
These sources were selected with the intention of identifying high-redshift analogs of the z < 0.8 
active galactic nuclei (AGN) with double-peaked [O III]A5007 lines, which might represent AGN out- 
flows or dual AGN. Lines of high-ionization potential are believed to originate in the inner, highly 
photoionized portion of the narrow line region (NLR), and we exploit this assumption to investi- 
gate the possible kinematic origins of the double-peaked lines. For comparison, we measure the 
[Ne V]A3426 and [Ne III]A3869 double peaks in low-redshift (z < 0.8) [O III]-selected sources. We find 
that [Ne V]A3426 and [Ne III] A3869 show a correlation between line-splitting and line- width similar to 
that of [O III]A5007 in other studies; and the velocity-splittings are correlated with the quasar Edding- 
ton ratio. These results suggest an outflow origin for at least a subset of the double-peaks, allowing 
us to study the high-ionization gas kinematics around quasars. However, we find that a non-neligible 
fraction of our sample show no evidence for an ionization stratification. For these sources, the outflow 
scenario is less compelling, leaving the dual AGN scenario as a viable possibility. Finally, we find 
that our sample shows an anti-correlation between the velocity-offset ratio and luminosity ratio of the 
components, which is a potential dynamical argument for the presence of dual AGN. Therefore, this 
study serves as a first attempt at extending the selection of candidate dual AGN to higher redshifts. 

Subject headings: galaxies: active - galaxies: nuclei - quasars: emission lines 



1. INTRODUCTION 

Mergers of gas-rich galaxies are likely to play a 
key role in the growth of supermassive black holes 
(SMBHs) through ac cretion, particularly in the trigger- 
ing o f quasar phases ([Sanders et al.l ll988: Trc ister et al.l 
12010) . Quasars may also represent an important phase 
in th e evolution of galaxies due to radiative feed - 
back (|Silk fc Reesl [1991 IKauffmann fc Haehnelti l200l . 
This sequence of events is a potential scenario link- 
ing the evolution of ga l axies and the g rowth of SMBHs 
(|Di Matteo et al.l 120051 : iHopkins et alJl2005[ ). Interest- 
ingly, both the merger and quasar phases can manifest 
themselves as active galactic nuclei (AGN) with double- 
peaked narrow emission lines in their spectra. Specif- 
ically, during a galaxy merger, when the SMBHs are 
separated by ~1 kpc and are act ively accreting as a 
dual AGN ([Comerford et al.ll2009bl ) they will each pro- 
duce Doppler-shifted emission from their narrow line re- 
gions (NLRs), resulting in a double-peaked profile in the 
integrated spectrum. In the feedback scenario, radia- 
tion from a si ngle quasar can dri v e bi-conical outflows 
of the NLR (lArribas et al.l 119961 : iVeilleux et all 120011 : 
iCrenshaw et al.ll2010afl. producing similar em ission line 
profiles (|Zheng et al.lll990t iFischer et al.ll20lH) . 

AGN whose spectra exhibit double-peaked narrow 
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emission lines (i.e. they are best fit by two compo- 
nents) represent a sub-class of the AGN population. Ex- 
tended, offset and double-peaked line profiles have been 
observed in the narrow emission lines of AGN since th e 
earliest studies of AGN NLRs ([Heckman et al.l 119811 ). 
These complex narrow line profiles are observed in both 
Type 1 and Type 2 AGN ([Crenshaw et alJl2010bf ). and 
they are generally most pronounced in forbidden transi- 
tion lines of high ionization potentials (I.P.s). In par- 
ticular, such observations have often focused on the 
5007A transition line of [O III] (I. P. = 35.15 eV) since 
it is a relatively intense emission line produced by the 
ionizing contin uum of AGN and is accessible i n opti- 
cal spectra (see I Veilleuxl 119911 and iWhittlel 119921 for ex- 
amples). However, there is even variation among the 
high-ionization lines, with those of the highest ioniza- 
tion potentials, such as [Ne III] (I. P. = 41.07 eV), and 
[Ne V] (L P = 97.16 eV), displaying the largest veloc- 
ity offsets dDe Robertis fc Osterbro"ck 1984: S turm et al.l 
[2002T:ISDOon fc Holtll2009T L~ 

In the dual AGN scenario, the two emission line peaks 
are produced by the orbital motion of two AGN within 
a single merger remnant galaxy. This interpretation is 
intriguing since SMBHs reside in the bulges of galax- 
ies, and dual SMBHs (kpc-scale separations) are a stage 
of galaxy mergers before the SMBHs coalesce. The 
existence of dual AGN has been confirmed observa- 
tionally in several serendipitous cases, most notably in 
ultra-lum inous infrared galaxie s as pairs of X-ray point 
sources (iKomossa et all l2003t iGuainazzi et al.l 120051: 
Hudson et al.l 120061: iBianchi et al.l 120081: iPiconcelli et al.l 
20101 : IKoss etaiTl2011l : IMazzarella et al.l 12012ft . SysT 



2 



Barrows et al. 



tematic searches for dual AGN in large spectro- 
scopic database s, such as the Sloa n Digital Sky 
Survey (SPSS: lAbazaiian et all f2009h and DEEP2 
([Newman et al.l 120121 ) have involved identifying AGN 



there is significant controversy over the role that galaxy 



with double emission l i ne systems dC omerfor d et al 
2009bl iWang et all [200l ILiu et alJ l2010bt ISmith et al 
2010t iGe et al.1 I2012T ). Promising results are being 



obtained through follow-up observati ons in the form 
of hig h-resolution optical imaging (|Comerford et all 
l2009aj). near-infrared (NIR) adaptive optics imag- 
ing (ILiu et al.1 l2010at iFu et al.l l2011at iRosario et al.1 
120111: IShen et all 1201 laHBarrows et all 120121) . spatially 
resolved spectroscopy (jMcGurk et al.l 120111: IFu et all 



20121). hard X-ray observations dComerford et all 120111 : 



Civano et all 1201 2t ILiu et alll2012al ). radio observations 
(|Fu et all 120 lib! ) . and diag nostics deduced statisti cally 
from longslit spectroscopy (jComerford et all 120121 ) . A 
sample of dual AGN will be useful in studying the con- 
nection between galaxy interactions /mergers and AGN 
activity (|Green et al.ll2010t ILiu et alll2011l l2012bf) . and 



for refining the gala x y merger ra t e (IConselice et al l2503l; 
Berrier et all 120061 : ILotz et all 120111 : iBerrier fe 
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In the outflow scenario, offset or double-peaked nar- 
row emission lines are often attributed to radially flowing 
NLR gas driven by energy from the AGN. Furthermore, 
this effect is also thought to produce a stratification of 
the NLR since the incident ionizing flux and electron 
density should diminish with increasing distance from 
the nuclear source. This will result in the production 
of different lines in varyi ng proportion s as a function 
of radius from the AGN (|Veilleuxiri99lD . Observation- 
ally, luminous quasars are known to have extended NLRs 
driven by energ y from the AGN cou pled to the interstel- 
lar medium gas (jBennert et al.ll2002l ) , and powerful radio 
galaxies often show complex, spatially extended NLRs 
with multiple components aligned along the radio axis. 
Models f or this alignment includ e reflection of the AGN 
emission (Tad hunter et allll988f) or material entrained i n 
a radio jet (jHolt et al II2003L l200l iKomossa et~alll2008l ). 
In AGN with sufficiently high Ed dington ratios, radiation 
pressure acting on gas and dust (lEverettll2007f) or a hot 
wind that entrains the NLR clouds ()Everett fc Murray! 
2007) are possible mechanisms capable of driving the 
line splitting. Whichever mechanism is the dominant 
driver of outflows in an individual source, they may 
represent cases of AGN feedback, which might be im- 
portant in quenching star formation in galaxies follow- 
ing a merger and in establishing the observed correla- 
tions between SMBH masses and host galaxy prop erties 
(jFerrarese fc Merrittl[2000l: iMarconi fc Huntll20q3l) . 

Whether AGN with double-peaked emission lines rep- 
resent dual AGN or powerful AGN outflows, they have 
proven useful for investigating several aspects of AGN 
evolution, including AGN triggering and feedback. Sys- 
tematic searches for double-peaked AGN have primar- 
ily utilized the [O III]A5007 emission line and have 
therefore been limited to below redshifts of z « 0.80 
since [O III]A5007 is not accessible in optical spectra 
at higher redshifts. However, there is evidence that 
at high redshifts galaxy mergers were more prevalent, 
and AGN outflows might have played an important role 
in the evolution o f galaxies at z ^> 1 , including mas- 
sive radio galaxies (Nesvad ba et al.1 120081 Furthermore, 
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2012j). Therefore, a sample of high red- 



shift double-peaked narrow line AGN will be important 
for investigating these questions. So far only one such 
candidate dual AGN has been identified above z ~ 0.8 
as a serendipitous discovery at z — 1.175 through double- 
peaked UV and optical emission lines, p articularly evi- 
dent in [Ne V]A3426 and [Ne III]A3869 (| B arrows et all 
I2012T ) . Motivated by this discovery, we have conducted 
a systematic search for additional AGN at high red- 
shift (z > 0.8) with double-peaked [Ne V]A3426 and 
[Ne III]A3869 emission lines. These emission lines have 
relatively high ionization potentials and trace the gas 
photoionized by the AGN continuum. Furthermore, they 
are likely to originate in the inner (more highly ionized) 
portion of the NLR, allowing for the identification of 
AGN driven outflows at 0.8 < z < 1.6. This study will 
also serve as a first attempt at extending the identifica- 
tion of candidate dual AGN to higher redshifts. 

In Section [5] we describe our parent sample and how 
we selected our final sample and low-redshift compari- 
son sample. In Section [3] we describe the systemic red- 
shifts we will use throughout the paper. In Section 0] 
we describe the general properties of our sample and 
our estimates of selection completeness. In Section [5] 
we investigate several correlations among the emission 
line properties that aid in illuminating the origin of the 
double-peaked line profiles. In Section |6] we describe the 
radio properties of our sample and compare them to our 
parent sample. In Section [7] we discuss the most likely 
physical mechanisms driving the line-splitting, particu- 
larly focusing on the scenarios of AGN outflows and dual 
AGN. In Section [8] we summarize our main conclusions. 
Throughout the paper we adopt the cosmological pa- 
rameters n A = 0.728, fl b = 0.0455, Q m h 2 = 0.1347, 
and Ho = 70.4 km s ?1 Mpc ?1 . This corresponds to 
the maximum likelihood cosmology from the combined 
WMAP +BAO+H0 results fro m the WMAP 7 data re- 
lease of IKomatsu et all (|2011D . 

2. GENERATING THE SAMPLE 
2.1. Parent Sample 

Our parent sample consists of archival spectra drawn 
from the quasar catalog of the SDSS Data Release 7 
(DR7) which is described in detail in ISchneider et all 
(2010). The typical resolution of the SDSS spectra is 
A/AA ~2000. In short, inclusion in the catalog requires 
luminosities brighter than Mj — —22.0, at least one 
emission line with FWHM > 1000 km s _1 or complex 
absorption features, apparent magnitudes fainter than 
i ss 15.0, and have highly reliable redshifts (see Section[3] 
for a discussion of the redshifts). We restricted the lower 
rcdshift limit of the sample to z > 0.80 to only include 
sources not in the parent samples of [O Hfl-selected 
doubl e-peaked emitters from the SPS S (IWang et all 
l2lMlLmet alll2010ri ISmith et al.llMTol: [Geet al.ll2012T) 
since our intention is to select sources which are not 
identifiable by their methods. Additionally, we required 
that at least [Ne V]A3426 be accessible in the SPSS 
wavelength range (3800 - 9200 A), which imposes an 
upper limit of z ~ 1.7. We did not make any selection 
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cuts based on the signal-to- noise ratios (S/N). This 
resulted in a parent sample of 39,876 sources. 

2.2. Initial Selection 

Our initial selection involved visually identify- 
ing quasars from the parent sample (Section 12. ip 
with detectable double emission line peaks in 
[Ne V]A3426 and/or [Ne III]A3869. These two lines 
were used since they are accessible in the SDSS optical 
wavelength range at z > 0.80, are relatively strong 
narrow lines in quasar spectra, and are not severely 
blended with any other strong lines. We did not 
require that two explicit peaks be detectable in both 
of those emission lines for two reasons: 1) the ionizing 
continuum may be such that [Ne V]A3426 is too weak 
to be detected whereas [Ne III]A3869 is detectable; and 
2) a difference in line ratios and/or velocity-splittings 
between the two lines may result in one pair being 
more blended than the other. Therefore, since the 
purpose of this analysis is to investigate the origin of the 
double-peaked emission lines, we did not want to exclude 
those sources which show variations among the line 
properties. While [O II]A3727 is another strong emission 
line accessible in most of our parent sample, we did not 
select sources based on this line since [O II]A3727 is a 
doublet (A3726,3729 A) with -200 km s" 1 separation 
between the transition wavelengths and is difficult to 
discern from true peaks, particula rly with the limi ted 
spectral resolution of the SDSS (see lSmith et al.|[2010l for 
a similar discussion). In sources where the double peak 
detections are ambiguous in each line, corresponding 
peaks in both [Ne V]A3426 and [Ne III]A3869 are needed 
for confirmation. Out of the parent sample, our vi- 
sual selection process resulted in 181 preliminary sources. 

2.3. Modeling the Spectra and Selection of the Final 

Sample 

To generate the final sample to be used in our subse- 
quent analyses, we modeled each spectrum yielded by the 
initial selection stage fSection l2.2|) in order to determine 
if a two-Gaussian model significantly improves the line 
profile fit over a single-Gaussian model. This modeling 
proceeded in two stages: 1) continuum modeling, and 2) 
emission line modeling. 

2.3.1. Continuum Modeling 

The continuum was modeled by masking all detectable 
emission lines and simultaneously fitting a power-law 
function (F\ — A~ QA ) for the underlying AGN con- 
tiuum radiation plus a pseudo-continuum of broad- 
ened Fe II emission lines from the empirical tem- 
plates of iTsuzuki et all (|2006l ) (A obs < 3500 A) and 
iVeron-Cettv et all ()2004f ) (A obs > 3500 A) which were 
developed from the spectrum of the narrow-line Seyfert 
1 galaxy I Zw 1. Due to the use of two separate Fe II tem- 
plates, the powerlaw function was allowed to have a 
break at 3500 A, and the Fe II normalization allowed 
to vary independently below and above the break wave- 
length. The Fe II pseudo-continuum is composed of 
many blended Fe II transitions which are generally be- 
lieved to originate in or near the classical broad line re- 



gion (BLR), and as such the redshift and broadening 
of the Fe II emission should be comparable to that 
of the broad emission lines. Therefore, in our model- 
ing the Fe II emission was fixed at the redshift of the 
quasar being modeled (see Section [3] for a detailed dis- 
cussion of the redshifts); we broadened the Fe II tem- 
plate by convolving with a Gaussian of FWHM~ conv 
where FWHM$ eII = FWHM? onv + FWHM? Zw x and 
FWHMi Zw i = 900 km s _1 . Use of this template 
necessarily limits the minimum FW HMpeii best-fit to 
900 km s _1 , though this was not a problem since these 
sources were in the quasar catalog based on the presence 
of broad lines with FWHM > 1000 km s _1 . The lower 
rest-wavelength end of the continuum+Fe II fitting win- 
dow was 2750 A and the upper rest-wavelength end was 
3950 A (just redward of [Ne III]A3869) if accessible, or 
otherwise the red end of the spectral coverage. This fit- 
ting window allowed proper detection of the Fe II emis- 
sion since it has a relatively large equivalent width near 
Mg IIA2800. We allowed FWHM FcU to vary in steps of 
100 km s , and in general, the solutions are in the range 
-2000-9000 km s _1 , though in most cases the quality of 
the fits were not strongly dependent on the broadening. 

2.3.2. Narrow Emission Line Modeling 

We fit each emission line ([Ne V]A3426 and 
[Ne III] A3869) with both a single and a double Gaus- 
sian model, and required that sources in our final sam- 
ple have fits to at least one of these two lines signifi- 
cantly improved by the double Gaussian model (A% 2 > 
Atiuof)- In seven individual line models, an ad- 
ditional broad component of line width FWHM w 
700 — 1500 km s _1 was required for a satisfactory 
fit. Three of the [Ne V]A3426 fits required such a 
component, where in two of those cases the broad 
component is consistent with the blue narrow peak 
(J074242.18+374402.0: FWHM = 1470 km s" 1 ; and 
J105035.57+190544.2: FWHM = 1220 km s^ 1 ), and in 
the third case it is consistent with the red narrow peak 
(J150243.93+281739.9: FWHM = 1320 km s^ 1 ). Four 
of the [Ne III]A3869 fits required a broad component, 
where in three of those cases it is consistent with the blue 
narrow peak (J105035. 57+190544.2: FWHM = 1520 
km s" 1 ; J105634.56+121023.5: FWHM = 820 km s -1 ; 
and J145659.27+503805.4: FWHM = 700 km s" 1 ) 
and in one those cases it is located between the blue 
and red peaks (J085205.91+183922.2: FWHM = 1020 
km s -1 ). We note that in only one of these sources 
(J105035. 57+190544. 2) is a broad component seen in 
both [Ne V]A3426 and [Ne III]A3869, and it is consistent 
with the blue peak. 

Properly accounting for the contribution of flux from 
Fe II emission is often important in determining the pro- 
files of the emission lines of interest in this wavelength 
range. In particular, there is a local peak in the broad- 
ened Fe II emission just blueward of [Ne V]A3426, be- 
tween 3390 A and 3410 A, which has the potential to 
complicate the isolation of the blue [Ne V]A3426 compo- 
nent. Therefore, careful attention was paid to the mod- 
eling in this region, and any sources for which the model 
was ambiguous were not admitted into the final sam- 
ple. However, there is no significant Fe II emission near 
[Ne III] A3869, and only a small local peak redward of 
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Fig. 1. — Four examples of quasar spectra yielded by our selection process shown in the quasar rest-frame (zsdss) wavelength range of 
3360 — 3950 A. The spectra have been smoothed by convolving with a Gaussian of a = 1 A. The flux densities, Fx> are m units of 10~ 17 
erg s — 1 cm -2 A -1 . In each of the four panels the power-law continuum has been subtracted from the spectrum, and the Fell template 
(gray, solid line), best-fit Gaussians (black, dotted lines), and the best-fit model sum (red, solid line) are shown. See Section [23] for details 
on fitting the models. These are examples for which double components are measurable in both [Ne V]A3426 and [Ne III]A3869. In all 
four cases we can only place upper limits on the presence of a second [O II]A3727 component (vertical dashed line). In each panel we 
have labeled the expected position of based on its rest-wavelength, however it is not always detected, and in no cases can we resolve 
double-peaks in H£. 
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Fig. 2. — Examples of two spectra from Figure [T] zoomed-in on 
[Ne V]A3426 (top) and [Ne III]A3869 (bottom). The spectra have 
been smoothed by convolving with a Gaussian of a = 2 A. In 
each panel the peak flux of the double Gaussian model has been 
normalized to unity, and the spectra are plotted in velocity-space 
with zero- velocity at the quasar redshift. The power-law continuum 
has been subtracted from the spectrum, and the model components 
are the same as described in Figure [l] 
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Fig. 3. — Distribution of reduced x 2 values for the best 
double-Gaussian fits for [Ne V]A3426 (black, solid line) and 
[NeIIl]A3869 (gray, dashed line). 

[O II]A3727 but with which it is not blended. 

From our spectral modeling, there are 38 sources in 
our sample for which we have measured robust double- 
peaks in both [Ne V]A3426 and [Ne III]A3869, 42 sources 
for which we have only robust [Ne V]A3426 measure- 
ments, and 51 sources for which we have only robust 
[Ne III] A3869 measurements. Thus, there are a to- 
tal of 131 double-peaked sources in our sample, mak- 
ing the fraction of double-peaked AGN detected in this 
manner ^0.3% of the parent sample. See Figure [I] 
for examples of the spectra in our final sample show- 



ing the best-fit models over the rest-wavelength range 
containing [Ne V]A3426 and [Ne III]A3869. Addition- 
ally, Figure [2] compares the individual [Ne V]A3426 and 
[Ne III]A3869 line profiles in velocity-space. The best- fit 
parameters of the emission line modeling are listed in Ta- 
ble Q] ([Ne V]A3426) and Table El ([Ne III]A3869). All er- 
rors correspond to la uncertainties and have been propa- 
gated throughout any further calculations. The reduced 
X 2 values for each double Gaussian model are listed in 
Tables Q] and [2] and shown as a distribution in Figure G3 
The emission line modeling was performed with SPEC- 
FIT (|KrissHl994f ). Our best-fitting Fe II FWHMs are 
listed in Table EH along with other relevant quasar prop- 
erties (redshifts, Mg IIA2800 FWHMs, and Eddington 
ratios) which will be described in subsequent sections. 

Throughout the rest of the paper, the velocity off- 
sets of the individual blue and red components will 
refer to the offsets from the systemic velocity. For 
[Ne V]A3426, these blue and red velocity offsets will be 
defined as AVp^v^biuc an d AV[ No y] !rod , respectively. For 
[Ne III]A3869, the blue and red velocity offsets will be 
defined as AVj NeII i] ,biuo and AVj NeIII ] >red , respectively. 
Blueshifts will correspond to positive velocities. The 
velocity-splittings will then be the velocity difference be- 
tween the lines, i.e. AV^oV = AVj NeV ],biue - AV[ Ne v],red 
and AVN e in = AVr NeIII ] ibluc - AV[ NoIII ] irod . The redshifts 
corresponding to the systemic velocity of each source are 
described in Section [3] 

2.4. Comparison Sample 

To generate a z < 0.8 comparison sample to be used 
in our analyses, w e examined the 89 [O III]-selected 
Type 1 sources in iSmith et al.1 (|2010l ) and measured 
the line properties of double-peaked [Ne V]A3426 and 
[Ne III] A3869 in cases for which two peaks are detectable 
in those lines. We chose this sample because, as opposed 
to other double-peaked [O III]A5007 samples, it includes 
Type 1 AGN which were selected from the SDSS quasar 
catalog. Therefore, it provides a z < 0.8 sample which 
can be used for direct comparison. For the measurements 
on the [O III]-selected sample, we selected those with ro- 
bust [Ne V] / [Ne III] double-peaks in the same way as for 
our high-redshift sample, and likewise modeled the spec- 
tra as described above. This resulted in 18 sources in 
the [O III]-selected comparison sample for which we have 
measured robust double-peaks in both [Ne V]A3426 and 
[Ne III]A3869, 4 sources for which we have only have 
robust [Ne V]A3426 measurements, and 20 sources for 
which we have only have robust [Ne III]A3869 measure- 
ments. Thus, there are a total of 42 double-peaked 
sources in our comparison sample. 

3. REDSHIFTS 

Some sections of our analysis require a knowledge of 
the individual quasar redshifts. For the sources in our 
sample, we have redshifts available from several differ- 
ent measurement techniques. In this section we describe 
those redshift estimates which are useful for our scientific 
interests. 

3.1. SDSS Redshifts 

Redshifts for spectroscopically-detected sources in the 
SDSS are determined by the 'spectrold' code described 
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TABLE 1 

Emission Line Properties for Double-peaked [Ne V]A3426 







Blue 








Red 








SDSS Name 


F A a 


AV b 


FWHM C 






AV b 


FWHM C 


Reduced 


X 2 (dof) 


J003159.87+063518.8 


5.4 + 4 


334 ± 124 


457 ± 268 


6 


9 + 4 


-128 + 6 


465 ± 203 


0.58 


(10) 


J014933. 86+143142.6 


4.7 + 2 


959 ± 55 


260 ± 97 


18 


7 + 4 


319 ± 57 


731 ± 209 


1.25 


(42) 


J015325. 74+145233.4 


28.3 ± 8 


504 ± 49 


403 ±119 


33 


9 + 11 


- 176 ± 68 


585 ± 249 


0.70 


(32) 


J021648. 36-092534.3 


7.0 + 2 


930 ± 35 


245 ± 50 


34 


9 + 5 


184 ± 41 


802 ± 131 


0.65 


(43) 


J021703. 10-091031.1 


16.4 + 4 


871 ± 35 


364 ± 86 


28 


+ 5 


155 ± 51 


671 ± 156 


0.82 


(30) 


J034222. 54-055727.9 


23.8 ± 6 


746 ± 73 


641 ± 206 


24 


4 + 6 


- 136 ± 83 


682 ± 176 


0.90 


(85) 


J035230.55-071102.3 


149.1 ± 100 


690 ± 270 


933 ± 267 


244 


1 ± 108 


- 15 ± 110 


788 ± 185 


0.72 


(79) 


J073408. 62+411901.1 


10.8 + 4 


589 ± 81 


364 ± 214 


20 


1 + 5 


-11+68 


469 ± 151 


1.95 


(35) 


J074242. 18+374402.0* 


9.9 ± 5 


198 ± 34 


236 ± 125 


10 


3 + 4 


- 116 + 30 


209 ± 74 


0.81 


(33) 


J074641. 70+352645.6 


38.9 ± 22 


1070 ± 327 


1306 ± 476 


22 


8 + 21 


124 ± 254 


961 ± 373 


0.84 


(58) 



Note. — All errors correspond to lcr uncertainties. * denotes that a broad component was necessary. This table is available in its 
entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 

a Flux densities, Fx, are in units of 10 — 17 erg s _1 cm -2 A -1 . 
b AV is velocity offset from the quasar redshift in km s — 1 . 
c FWHM are the observed line widths in units of km s — . 



TABLE 2 

Emission Line Properties for Double-peaked [Ne III]A3869 







Blue 






Red 








SDSS Name 




AV b 


FWHM C 


p A a 


AV b 


FWHM C 


Reduced 


X* (dof) 


J000531. 41+001455.9 


28.0 + 7 


146 ± 53 


372 ± 125 


11.2 + 5 


- 179 ± 32 


125 ± 49 


2.15 


(30) 


J003159.87+063518.8 


6.2 + 11 


446 ± 139 


223 ± 278 


18.1 ± 12 


61 ± 232 


560 ± 274 


2.07 


(31) 


J004305.92-004637.6 


7.8 + 6 


395 ± 74 


236 ± 276 


9.5 + 6 


111 ± 55 


209 ± 151 


1.20 


(32) 


J004312. 70+005605.0 


9.3 + 2 


215 + 31 


215 + 65 


16.6 + 2 


- 140 ± 13 


212 ± 24 


1.32 


(33) 


J014822.62+132142.7 


37.8 + 15 


883 ± 78 


1048 ± 438 


27.6 + 6 


- 376 ± 102 


779 ± 149 


0.83 


(47) 


J015734.24+003405.5 


17.1 + 5 


806 ± 50 


260 ± 93 


70.6 + 12 


-69 + 42 


672 ± 131 


1.05 


(37) 


J021648.36-092534.3 


8.9 + 2 


538 ± 27 


185 ± 50 


35.8 + 5 


-4 + 36 


546 ± 95 


1.06 


(24) 


J023234.33-091053.0 


51.4 + 17 


1462 ± 157 


1565 ± 506 


34.4 ± 10 


31 ± 151 


1121 ± 238 


0.73 


(66) 


J034222. 54-055727.9 


28.8 + 9 


480 + 


925 ± 302 


58.4 + 7 


- 323 ± 37 


656 ± 81 


1.09 


(78) 


J035230.55-071102.3 


120.0 + 42 


382 ± 68 


433 ± 135 


201.2 + 40 


-70 + 33 


379 ± 57 


0.95 


(32) 



Note. — All labels are the same as in Table [T] This table is available in its entirety in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding its form and content. 



TABLE 3 

General Quasar Properties for Our Sample 



SDSS Name 




ZSDSS 






ZMgll 




FWHM Fen a 


FWHM M Gii a ' b 


fEdd C 


J000531 


41+001455.9 





9918 ± 


0009 





9931 ± 


0018 


9000 


5424.76 + 1192.59 


0.109 


J003159 


87+063518.8 


1 


0921 ± 


0010 


1 


0935 ± 


0016 


1000 


2475.19 + 1745.19 


0.085 


J004305 


92-004637.6 





8482 ± 


0015 





8488 ± 


0015 


8400 


2828.74 + 247.77 


0.094 


J004312 


70+005605.0 





9036 ± 


0010 





9047 ± 


0014 


8500 


6223.22 ± 283.39 


0.105 


J014822 


62+132142.7 





8767 ± 


0018 





8759 ± 


0019 


8000 


6063.22 ± 134.05 


0.146 


J014933 


86+143142.6 





9024 ± 


0013 





9027 ± 


0013 


4500 


2969.71 ± 244.58 


0.172 


J015325 


74+145233.4 


1 


1755 ± 


0020 


1 


1762 ± 


0021 


5100 


5219.56 + 553.29 


0.202 


J015734 


24+003405.5 


1 


0834 ± 


0013 


1 


0825 ± 


0020 


3300 


6846.37 + 1636.32 


0.109 


J021648 


36-092534.3 





8834 ± 


0013 





8834 ± 


0014 


2700 


4288.12 + 632.18 


0.119 


J021703 


10-091031.1 





8752 ± 


0013 





8757 ± 


0015 


2700 


2598.31 + 432.35 


0.289 


Note. 


- This table 


is available 


n its entirety in a 


machine-readable form in the online journal. A 



portion is shown here for guidance regarding its form and content. 

a FWHMs are in units of km s — 1 . 

b Values for FWHM MgII are from the SDSS DR7 Catalog of Quasar Properties llShen et al.H2011bl) . 
c Eddington ratio calculations are described in Section 1-5.21 
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Fig. 4. — Distribution of (zsdss - z Mg ii)/(1 + zsdss) for the 
parent sample (top) and our sample (bottom). The mean and 
sigma are shown for each distribution. Note that the distributions 
are similar for both samples. 



in iStoughton et afl (|2002l ). In short, two separate red- 
shifts are determined, an emission line rcdshift (zel) and 
a cross-correlation redshift (zxc), and the final value 
adopted by the code is the redshift solution with the 
highest confidence level (though each individual value 
is stored and available through the archive). Addi- 
tionally, a small fraction of the quasar redshifts were 
re-determined after v isual inspection as described in 
ISchneider et~ai1 ((2010). zel is measured by the identifi- 
cation of common galaxy and quasar emission lines with 
known rest- wavelength values, and therefore for our sam- 
ple the highest confidence zel values are almost exclu- 
sively determined by sets of the strong quasar emission 
lines C IVA1549, C III]A1909, Mg IIA2800, [O II]A3727, 
H<5, H/y, and H/3. Likewise, the template which pro- 
vided the highest confidence zxc val ues for our sample is 
the co mposite quasar template from IVanden Berk et al] 
(|2001[ ). For 107 of our sources, zxc and zel are consis- 
tent with each other within their errors. For the remain- 
ing 24 redshifts, 19 are from zxci 4 are from Zel, and 1 
was determined by hand. We note that for those sources 
with disagreeing values of zxc and zel, the poorer of the 
two values is clearly incorrect and not usable. Through- 
out the rest of this paper, we refer to the final redshifts 
produced by the 'spectrold' code as zsdss, and these 
values are listed in Table [3] 

3.2. Redshifts from Mg 77A2800 

In principle, individual emission lines can provide in- 
dependent redshift estimates corresponding to the phys- 
ical regions where those lines originated, e.g. the broad 
emission lines provide redshifts for the BLR. For our pur- 
poses, it is useful to know the redshift of the central 
SMBH (under the assumption of a single, active SMBH 
in the host galaxy) which can potentially be traced by the 
BLR if the gas is virialized. Fortuitously, given the red- 
shift range of our sample, for all of our sources we have 




Fig. 5. — Top: Distributi on of the parent SDSS quasar 
sample (Schneider et al. 2010) in the relevant redshift range. 
Bottom: Distribution of redshifts for our final double-peaked 
[Ne V]/[Ne III] sample developed as described in Section [2] The 
cut-off below z = 0.80 in our sample is ent irely the result of our 
initial cut in redshift space (see Section [2TJ , and the highest red- 
shift for sources found through our selection process is z = 1.53. 
The redshifts are described in Section [3] 



spectral access to the broad Mg IIA2800 emission line, 
which has been sho wn to be virialized in th e potential of 
the central SMBH (jMcLure fc Jarvisll2002D . Therefore, 
we have used the Mg IIA2800 emission line models from 
the ca talog of SDSS DR7 quasar properties (tShen et al.1 
1201 lbf ) to obtain estimates of the BLR redshifts (zMgii)- 
However, we would like to note several potential techni- 
cal difficulties involved in measuring the Mg IIA2800 line 
centroid which limit its application to our work. These 
difficulties include blending with Fe II emission, UV 
absorption features, and the often ambiguous presence 
of a narrow emission line component (which is actu- 
ally the doublet Mg IIAA2797, 2802). We also note 
that the Mg IIA2800 line centroids may be compli- 
cated if the low-ionization broad emission lines have 
the characteristic double-peaked profile of disk e mitters 
(jEracleous fc Halpernl l200l iStrateva et all 120031 ). Fur- 
thermore, since we are considering the scenario of two 
AGN within the same host galaxy, we must consider the 
possibility that the broad Mg IIA2800 line is a blend- 
ing of two BLRs from two Type 1 AGN. Though the 
components of a dual AGN would not be close enough 
for the line-splitting to exceed the FWHM (several thou- 
sand km s -1 ) and produ ce explicit double-peaked pro- 
files (jShen fc Loebll2010l ). two broad Mg IIA2800 com- 
ponents separated by several hundred km s _1 will still 
result in a relatively broadened, and possibly asymmet- 
ric, Mg IIA2800 profile, thereby complicating the cen- 
troid measurement. We have visually inspected the 
Mg IIA2800 line profile for each source in our sample 
in order to characterize their structure. There are sev- 
eral sources which show evidence for asymmetric struc- 
ture, though it is generally difficult to discern if any of 
the profiles contain two broad components as would be 
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the case if there are dual BLRs. Though many of those 
redshifts have rather large errors, a subset of them have 
robust line centroids with percent errors less than 1% (73 
sources) and percent errors less than 0.1% (52 sources). 
The Mg IIA2800 redshifts are listed in Table [3] 

3.3. Comparison of zsdss and ZMgii 

Figure0]shows the distribution of (zsdss— ZMgii)/0-+ 
zsdss) where the Mg IIA2800 redshifts appear to be sys- 
tematically larger than the cross correlation redshifts in 
both our sample and the parent sample. This same ef- 
fect is clearly a p paren t in the SDSS redshift analysis of 
iHewett fe Wildl ([2010), which suggests that we are see- 
ing the same systema tic trend in our analysis. While 
IHewett fc Wildl (|2010l ) provide redshifts which correct 
for this systematic offset in a statistical sense, these 
corrections may not be appropriate for some individual 
quasars. Therefore, we do not use these improved red- 
shifts since our sample only contains 131 sources. 

Based upon our comparison between zsdss and ZMgii, 
for the rest of our analysis we have chosen to use the 
zsdss values and their associated errors. These SDSS 
redshifts will serve as the systemic velocities for our 
sources. All quoted line-splittings between the two line 
peaks (AV[N e v] an d AV[ NeIII ]) and velocity offsets for the 
blue (AV [NeV]jbluc and AV [NcIII ] ibluc ) and red (AV[ Ne v],red 
and AVfNeinLred) components are relative to these SDSS 
redshifts. The range of redshifts in our sample is z = 0.80 
(low-limit cutoff discussed in Section I2.1[) to z = 1.53 
(highest redshift source in our sample). In Section IS~4l 
for which our analysis is heavily dependent upon the 
choice of redshift, we will also us e th e robust ZMgii mea- 
surements mentioned in Section 13.21 for comparison. 

4. GENERAL PROPERTIES OF THE SAMPLE 

In this section we describe several of the general prop- 
erties of our sample. We compare the distributions of 
these properties (redshifts, S/N, and velocity offsets) to 
the parent sample or our low-redshift comparison sam- 
ple in order to show our selection biases. Additionally, 
we discuss the results of our completeness estimates with 
respect to several of these and other properties. 

4.1. Distributions 

The SDSS redshifts we have adopted for our sample 
(Section are shown and compared to the parent sam- 
ple in Figure [5] The redshift distribution shows a peak 
near the low redshift cutoff of the sample, with a gradual 
decline out to higher redshifts, in contrast to the increas- 
ing population of the parent SDSS quasar sample. This 
strong dependence on redshift may reflect a dependence 
on equivalent width or a luminosity bias, which we dis- 
cuss further below. 

Figure |5] shows the distribution of the continuum S/N 
per pixel, where the mean values are 18.5 and 14.5 in 
the wavelength regions adjacent to [Ne V]A3426 and 
[Ne III]A3869, respectively. Compared to the parent sam- 
ple, our sample shows a deficit at the low S/N regime, 
and a stronger tail out to higher S/N values. These 
distributions show that our selection is biased toward 
spectra of good quality which reflect the fact that de- 
tection of double-peaked structure in [Ne V]A3426 and 
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Fig. 6. — Distribution of continuum S/N per pixel near the 
[Ne V]A3426 (black, solid line) and [Ne III]A3869 (gray, dashed 
line) emission lines for sources in our sample. For each distribution, 
the S/N values were calculated over a set of featureless wavelength 
regions near the respective emission lines. 

[Ne III]A3869 is sensitive to the S/N and generally diffi- 
cult at S/N<5. Furthermore, the larger mean value of 
the continuum S/N near [Ne V]A3426 is reflective of the 
fact that [Ne V]A3426 is general ly a weaker line than 
[Ne III1A3869 in AGN spectra (|Osterbrock fc Ferlandl 
2006) and therefore a stronger signal is required for the 
detection of double-peaked structure. The dependence 
of our selection completeness on S/N will be addressed 
below. 

Since this study is motivated by the identification 
of high-redshift analogs of AGN with double-peaked 
[O III]A5007 lines, we compare the distribution of 
line properties of our sample to those of the [O III]- 
selected comparison sample. The SDSS spectral res- 
olution, S/N, and the intensity of the emission lines 
determine whether or not we are able to reliably de- 
tect double-peaks at a given separation in velocity- 
space. [Ne V]A3426 and [Ne III]A3869 are fainter 
th an [O III1A5007 in AGN sp ectra by a factor of ~ 
4 (|Ferland fc O sterbrock 1986), which means that se- 
lection based on [Ne V]/[Ne III] will be biased to- 
ward larger line-splittings relative to [O III]A5007. 
For example, double-peaked samples selected through 
[O III]A5007 have sources with line splittings down to 
-200 km s" 1 , and the mean of th e [O III]A5007 line - 
splittings in the Type 1 sample of S mith et al.1 (|2010f ) 
is —420 km s _1 . On the other hand, the mean of the 
[Ne V]A3426 and [Ne III]A3869 line-splittings of our sam- 
ple is ~700 and —700 km s _1 , respectively. 

Figure [7] shows the distribution of velocity off- 
sets from the systemic redshift for [Ne V]A3426 and 
[Ne III]A3869 in both our sample and the [O III]-selected 
sample. While the measurements on the comparison 
sample were performed in the same way as for our sample 
(Section [2]), these sources were selected from a sample of 
previously determined double-peaked AGN. Therefore, 



Outflows and Candidate Dual AGN at z = 0.8 -1.6 



1.0 

0.8 

«N 0.6 

0.4 

0.2 
0.0 
1.0 

0.8 

«N 0.6 

0.4 

0.2 
0.0 



[NeV]: This Sample 
AV blue : +740 km s"' + 




[NeV]: from S10 

AV blue : +320 km s"' \ 
AV red : -120 km s" 1 + 




[Nelll]: This Sample 

AV blue : +470 km s' 1 
AV red : -120 km s" 




[Nelll]: from S10 
AV blue : +300 km s" 1 



190 km s" 



-1000 1000 2000 

Velocity Offset (km s" 1 ) 



-1000 1000 2000 

Velocity Offset (km s -1 ) 



Fig. 7. — Distribution of velocity offsets for double-peaked AGN with positive velocities indicating blucshifts. 'Blue' component offsets 
(AVfNeVLblue an< ^ ^^[Nelll] ,blue) are shown as horizontally-hatched blue lines, and 'red' component offsets (AVrN e v],red an< ^ AV^Nenii red) 
are shown as diagonally-hatched red lines. The peak value in each panel has been normalized to unity. Top: Sources selected through 
[Ne V]A3426 and [Ne III] A3869 emission line profiles. Bottom: Distribution of velocity-offsets for the [O IIl]-selected Type 1 AGN from 
the lSmith et al.l H2010T ) (S10) sample. All distributions are shown separately for [Ne V]A3426 (left) and [Ne III]A3869 (right). The average 
velocity-offsets are shown for each distribution. The black solid vertical line in each panel represents zero velocity offset. 



we were able to more confidently identify two peaks in 
[Ne V]A3426 and [Ne III]A3869 at relatively smaller AVs 
for the [O III] -selected sample, as shown in Figure [7J A 
Kolmogorov-Smirnov (KS) test indicates that the blue 
components of our sample come from a different distri- 
bution than those of the [O III]A5007-selected sample 
(D[NeV] — 0.56, -DrNeiii] = 0.43, and probabilities that 
they come from the same distribution P[n c v] — 3 x 10~ 5 , 



[NcIII] 



7 x 10 ). A KS test also reveals a (less signif- 



icant) difference in the red components (-D[n c v] = 0.43, 
-D[Neiii] = 0.33, and probabilities of P[n c v] = 2.3 x 10~ 3 , 



[NcIII] 



5.6 x 10" 



This indicates that we are find- 



ing similar double-peaked narrow line sources to those 
selected through [O III]A5007 but that we are prefer- 
entially missing a portion of those with relatively small 
AV'b. 

4.2. Selection Completeness 

To estimate what fraction of sources with double- 
peaked [Ne V]A3426 and/or [Ne III]A3869 we are miss- 
ing due to the potential selection biases discussed 
above, we have estimated our selection completeness 
through simulations. Specifically, we generated arti- 
ficial spectra designed to mimic SDSS spectra in the 
[Ne V]/[Ne II I] wavelength ra nges (a similar test was 
performed by ILiu et all I2010bl for their [O III]-selected 
sample). Using the SDSS spectral dispersion, we mod- 
eled the continuum as a power-law and the emission lines 
as Gaussians with randomly assigned values for the pa- 
rameters of equivalent width (EW), FWHM, continuum 
S/N per pixel, and line offsets from the systemic velocity. 
The ranges of allowed values for those parameters were 



made uniform and at least slightly larger than the dis- 
tributions of our final sample in order to ensure that we 
sampled the relevant parameter space for the analysis. 
While the EW distribution of our sample only extended 
to ^35 A for the blue and red components, we allowed 
our simulations to have EWs of up to 50 A for the indi- 
vidual lines (100 A total EW). The simulated FWHMs 
ranged from 100 to 1800 km s _1 , ^200 km s _1 larger 
than the largest FWHM of our sample. The simulated 
continuum S/N values ranged from 2 to 33, just larger 
than the range of our sample (see Figure [5]) . Finally, 
the simulated line offsets ranged from to 2000 km s -1 , 
^500 km s _1 larger than the largest line offsets of our 
sample. We generated a total of 10,000 spectra, with 
2,000 of them having double-peaked emission lines, and 
the other 8,000 having single-peaked emission lines. Out 
of these simulated spectra, randomly distributed, we se- 
lected a final sample of double-peaked emission lines in 
the same manner as described in Section [5J 

We see a positive trend between completeness and EW, 
and find that we are highly incomplete over the EW 
range of our sample. For example, we estimate a com- 
pleteness of just ~1.2% for a total EW of 10 A and 
~4.5% for a total EW of 30 A. We see a negative trend 
between completeness and FWHM which can be inter- 
preted as the result of peak blending at large line widths. 
We see a moderate positive trend between completeness 
and S/N, with a completeness of ^42% at S/N=5 and 
^50% at S/N=30. Finally, we see a strong dependence 
on velocity-splitting, with a completeness of ~55% at 
AV = 9 00 km s" 1 and ~3 4% at AV = 500 km s _1 . In 
contrast. ILiu et all (|2010bD find a completeness of ~75% 
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at AV = 900 km s" 1 and -50% at AV = 500 km s" 1 
for their [O III]A5007 sample. 

Over co mparable paramete r ranges, we are less com- 
plete than iLiu et all (|2010b| ) for all of the parameters 
(~50-80% of their completeness), with the exception of 
FWHM, for which we have an equivalent or slightly 
greater completeness at the large FWHM end of their 
distribution (—900 km s^ 1 ). This is likely due to the fact 
that we allowed for a larger range of line offsets result- 
ing in our ability to identify broader double peaks. The 
s maller completenes s in our selection compared to those 
of lLiu et al.l (|20 10b) is likely due to the relatively smaller 



equivalent widths of [Ne V] / [Ne III] allowed in our simu- 
lations which is reflective of the weaker intensities com- 
pared to [O III]A5007. Of the parameters investigated, 
we find that our completeness is most s trongly dependent 
on and most drastically different from ILiu et al.l (|2010b[ ) 
in the velocity-offsets (Figure[7|). This result is again con- 
sistent with the notion that we are preferentially missing 
a portion of those with relatively small AVs. 

Since false double-peaks are generated by noise, false 
positives are most likely to occur in low S/N spectra 
in which random noise peaks are difficult to discern 
from the emission line signal. However, our conserva- 
tive selection criteria effectively required that the S/N 
be sufficiently large, as can be seen in Figure with 
S/Npss e Y] < 5 for only 4 sources (3%) and S'/iV[ No ni] < 5 
for only 11 sources (8%). This is reflected in the fact that 
we did not select any false positives in our completeness 
analysis. 

5. CORRELATIONS AMONG SAMPLE 
PROPERTIES 

In the following section, we examine several trends re- 
lated to the velocity offsets of the individual blue and 
red components from the systemic redshift and line- 
splittings between the blue and red components. We 
are interested in uncovering the prevailing mechanism(s) 
producing the double-peaked emission lines in our sam- 
ple. In particular, we investigate the following four rela- 
tionships: line-splitting vs line width, line-splitting vs 
quasar Eddington ratio, [Ne V]A3426 velocity offsets 
vs [Ne III]A3869 velocity offsets, and blue/red velocity- 
offset ratio vs blue/red luminosity ratio. In the following 
analyses we combine the [Ne V]/[Ne III] measurements 
of our sample with those of the [O III]-selected sample in 
order to increase both the sample size and the dynamic 
range of velocity offsets and line-splittings. 

5.1. A Correlation Between Line-Splittings and Line 

Widths 

In AGN with emission lines offset from the systemic 
velocity, there is often an observed positive correlation 
between the line peak offset (blueshifts being positive 
offsets) and line width. For example, in samples of 
so-called 'blue outliers' (AGN with [O III]A5007 lines 
blucshifted by > 100 km s _1 ), there is an observed pos- 
itive correlation between the blueshi ft and line width 
of [O III]A5007 dKomossa et all 12008ft . A similar trend 
was also seen by ILiu et al.l (|2010bl ) for their sample of 
Type 2 double-peaked [O III] A5007 sources, with the 
offsets being the line-splitting between the two peaks. 
In a similar fashion, we have plotted the line splittings 
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Fig. 8. — Plot of line width (FWHM) vs. line-splitting for 
[Ne V]A3426 (top panel) and [Ne III]A3869 (bottom panel). In each 
panel, the 'blue' and 'red' systems are labeled with blue circles and 
red triangles, respectively. Filled data points correspond to our 
[Ne V]/[Ne IIl]-selected sample, while unfilled points correspond to 
the [O IIl]-selected sample. The Spearman rank coefficients (rg) 
and probabilities of a null hypothesis (P) are listed for the 'blue' 
and 'red' components. The dashed line is the one-to-one relation. 
The average X and Y errors are indicated in the lower right cor- 
ner of each panel. Our completeness tests indicate that we are 
increasingly incomplete toward large-FWHM/small-AV (i.e. the 
upper left corners), but that we are increasingly complete toward 
small- FWHM/large- AV (i.e. the lower right corners). 



against the line widths (FWHM) for the double-peaked 
[Ne V]A3426 and [Ne III]A3869 lines in both our sam- 
ple and in the [O III]-selected sample (Figure [H]). For 
both [Ne V]A3426 and [Ne III]A3869, statistically signifi- 
cant positive correlations between the line-splitting and 
FWHM are evident. The correlations are strong and the 
line-splittings are generally larger than the li ne widths, 
and bo th of these results were also found by ILiu et all 
(201QU with respect to [O III]A5007 emission lines. 

The possible physical interpretations of these trends 
often includes an ionization stratification of the NLR 
in which the higher ionization lines originate closer 
to the ionizing radia tion source (jZamanov et al.l 120021 : 
iKomossa et al.ll2008f ). In such a scenario, if there is a 
radially decelerating outflowing component of the NLR, 
then the higher ionization lines are accelerated to higher 
velocities. The emission lines produced in the inner por- 
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Fig. 9. — Plot of Eddington ratio versus velocity-splittin g (to p), blue component offsets (middle) and red component offsets (bottom) for 
[Ne V]A3426 (left) and [Ne III]A3869 (right). See Section [572] for details on calculations of /Edd- To be consistent among each subsample 
while accommodating the negative velocity-offsets of red components, the individual velocities for each distribution have been shifted by 



a constant, a, such that the mean of that particular distribution is equivalent to 1000 km s 



a = 1000 - (AV). In each panel the 



horizontal gray, dashed line represents AV = km s — 1 for that particular distribution. Sources at 0.8 < z < 1.6 ([Ne V]/[Ne IIl]-selected) 
are plotted as filled triangles, and sources at z < 0.8 ([O IIl]-selected) are plotted as unfilled triangles. The Spearman rank coefficients 
and probabilities of the null hypothesis are given for each of the sub-samples and for the combined samples. The red solid lines in each 
distribution are the best linear fits weighted by the X and Y errors (obtained with the 'FITEXY' routine) for the combined sample. 



tions of the NLR will also be more broadened as their 
motion w ould be dominated by th e bulge gravitational 
potential (jNelson k Whittle! [1996). Since the ionizing 
potentials of [Ne V]A3426 and [Ne III]A3869 {LP. = 
41.07 and 97.16 eV, respectively) are larger than that 
of [O III]A5007 {LP. = 35.15 eV), we might expect 
the correlation between line-splitting and FWHM to be 
stronger among these emission lines. This is also ex- 
pected based on the higher critical electron densities 
for collisional de-excitation (n cr ;t) of [Ne V]A3426 and 
[Ne IH]A3869 (n crit = 1.3xl0 7 cm- 3 andn m( =9.5xl0 6 
cm~ 3 , respectively) compared to [O III]A5007 (rt C rit = 
6.8 x 10 5 cm~ 3 ) since electron densities will increase 
toward the nuclear region. Indeed, we find that the 
strengths of these correlations are generally comparable 
to or stronger than tho se for [O III]A5007 in the sam- 
ple of iLiu et~al"1 (|2010bt ) and in the narrow line Seyfert 



1 (NLS1) sample from iKomossa et~ai] (|2008[ ). These 
trends in our sample might be more reflective of the even 
stronger trends amon g the subset of 'blue-outliers' from 
IKomossa et"aL1 (|2008[ ). 

5.2. A Correlation Between Velocity-Splitting and 
Eddington Ratio 

Motivated by the outflow interpretation often used for 
samples of double-peaked emission line AGN, we investi- 
gate the relationship between the velocity-splittings and 
the quasar Eddington ratios (/Edd = ^boi/^Edd)- We 
calculated ^Edd using the standard derivation of the Ed- 
dington luminosity, LEdd = 47rcG-MBHMe/ (T T, where G 
is the gravitational constant, fi e is the mass per unit 
electron, and or is the Thomson scattering cross-section 
(|Krolildll999D . M BH was estimated for each Type 1 AGN 
using the SMBH mass-scaling relationships from broad 
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Fig. 10. — Top: Plots of AVb lue ,[ NeV ] versus AVbi ue ,[Nelll] (A) and AV^-edJNeV] versus AV red [ NeIII ] (B). The dashed line is the one-to-one 
relation. Each velocity has been shifted by the constant a where a = 1000 — (AVjjv e y]). In each panel the horizontal gray, dashed line 
represents AVj NeV ] = km s — 1 and the gray, shaded region represents the average of our sample. [Ne V]/[Ne IIl]-selected sources are 
plotted as filled circles while sources selected via [O III]A5007 arc plotted as open circles. Additional sources are plotted as stars: extended 
NLRs (SDSS J110851. 03+065901.4 and SDSS J135646.10+102609.0) and an unresolved NLR (SDSS J105052. 46+083934. 71 fromlFu et al. 
(2012) are indicated by solid green a nd purple stars, respectively, [OlII]-selected dual AGN SDSS J171544. 05+600835. 7 (Comcrford et al. 
120111) and SPSS J15 0243.1+111557 tFu et al. 2011b) are indicated by blue stars, and a [Ne V]/[Ne IIl]-selected dual AGN candidate from 
IBarrows et al.l H2012T ) is indicated by a red star (CXOXBJ142607. 6+353351). Bottom: Histograms of AVrNeV],Mue — AVjuelll] ,blue (C) 
and AVj N8 v], re d - AVj NcII i], ro d (D). The dashed line in the bottom panels represents AV [Nc v] ib lu C = AVj Ne iii],blue (C) and AV[ Ne v],red = 
AV[Nelll],red (D). The mean values are shown for each distribution. 



emission line widths and monochromatic luminosities: 
FWHMylp and L5100 A (z < 0.8) or FWHM MgU and 
L3000 A (z > 0.8) from McClure and Dunlop (2004). The 
line widths are m ostly from t he SD SS DR7 catalog of 
quasar properties iShen et all (|2011b[ ). though we visu- 
ally inspected and re-fit several of them for which we were 
able to improve on the models. For each AGN with an 
estimate of Mbh, we estimated Lboi from the monochro- 
matic luminosities L5100 A (z < 0.8) or L3000 A (z > 0.8) 
and the bolometric corrections of Richards et al. (2006). 

Figure [9] shows f^dd plotted against three quantities: 
velocity-splittings, blue velocity-offsets and red velocity- 
offsets. For both [Ne V]A3426 and [Ne III]A3869, a Spear- 
man rank test reveals statistically significant correlations 
between f^dd and each of these quantities. When com- 
paring [Ne V]A3426 and [Ne III]A3869, it is apparent that 
the correlations are stronger for [Ne V]A3426 than for 
[Ne III]A3869 in all correlations tested. Additionally, for 
both [Ne V]A3426 and [Ne III]A3869 the correlations are 
strongest between f^dd and the blue velocity-offsets. We 



note that the low velocity-offsets are dominated by the 
[O III]-selected sources, which simply reflects the selec- 
tion bias discussed in Section 21 If the velocity-splittings 
represent a NLR outflow in a significant number of these 
sources, then this result suggests that the radiation gen- 
erated by the quasar accretion rate may play a crucial 
role in driving the line-splittings. Additionally, it ap- 
pears that the blue component is more strongly domi- 
nated by outflowing material than the red component, 
which has implications for obscuration. 

5.3. Ionization Stratification 

The difference in LP. between [Ne V]A3426 and 
[Ne III]A3869 is ^56 eV and can be used to provide 
insight on the gas dynamics in different regions of the 
NLR. Since the line offsets have been shown to cor- 
relate with the quasar Eddington ratios (Section l5.2p . 
suggesting an outflowing component, we would like to 
search for evidence of an ionization stratification. Fig- 
ures [TOjA and [TUB show AV[NeV],blue verSUS AV[NoIII],blue 
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Fig. 11. — Plot of the velocity-offset ratio (| AVjl/ AV,-|) versus the luminosity ratio (L b /L r ) shown separately for [Ne V]A3426 (left) 
and [Ne III]A3869 (right). The solid black line is linear best fit weighted by the X and Y errors found using the 'FITEXY' routine. The 
solid red line is the theoretical relation assuming equivalent Eddington ratios for the blue and r ed co mponents (ei,2 = 1); and the dashed 
black line is the theoretical relation for a bi-polar outflow, both of which are discussed in S ecti on | 5. 41 The colored stars represent the same 
extended/unresolved NLRs, confirmed dual AGN and candidate dual AGN as in Figure ITul and have the same symbols. In each panel 
the Spearman rank correlations coefficients and probabilities of no correlation are shown for our sample, this comparison sample, and the 
combined ('All') sample. 



and AV[N C v],rcd versus AV[N e ni],red for all sources in our 
sample and the [O III]-selected sample for which there 
are detectable double-peaks in both lines. 

Though the true nature of the gas kinematics in NLRs 
is likely to be complex, under the simplest pictures of 
NLR ionization stratifications it is generally not expected 
that lines of lower ionization potential will show blue 
velocity-offsets larger than those of greater ionization po- 
tentials. Indeed, from Figure ITUR it can be seen that 
all of the points are consistent with having equivalent 
blue velocity-offsets in the two lines, or otherwise with a 
larger blue offset in [Ne V]A3426. This is also apparent 
in Figure [Top, with a mean AVj Ne v],biue ~ ^ ^Nelil] .blue 
value of 195 km s _1 . Accounting for the la uncertainties, 
from our sample 15 sources (40%) are consistent with 
equivalent blue velocity-offsets for both [Ne V]A3426 and 
[Ne III]A3869, and 23 sources (60%) have significantly 
larger [Ne V]A3426 blue velocity offsets, consistent with 
the presence of an ionization stratification. Compara- 
ble fractions are seen in the [O III]-selected sample. 
The sources in the small velocity-offset portion of the 
plot are dominated by the [O III]-selected sample, which 
again reflects the selection bias discussed in Section fJJ 
In contrast, the difference in the red offsets has a nar- 
rower distribution and a mean value more consistent with 
AV[NcV],rcd = AVjNeurj.red (Figures [TDB and [TUP). This 
suggests that the red systems are generally less stratified 



and have a different origin. 

Also plotted are several additional double-peaked 
sources for which follow-up observations have been 
obtai ned, including extended and u nresolved NLRs 
from IFu et all <f2012h dual A GN from IFu et all (12012ft 
and iComerford et al.1 (12011ft. a nd a candidate dual 
AGN from (jBarrows et all I2012f >. Note that the ex- 
tended/unresolved NLRs and the dual AGN are [O III]- 
selected, and the candidate dual AGN was selected based 
on [Ne V]A3426 and [Ne III]A3869. With the exception 
of the dual AGN SDSS J150243.1+111557, these addi- 
tional sources are all consistent with no apparent ioniza- 
tion stratification in either the blue or red system by the 
la criteria we have used above. 

5.4. Dynamical Relation 

In this section we examine a d y namic al argument, orig- 
inally proposed by IWang et al.l (|2009|) . for the presence 
of dual AGN in our sample. Under the simplest pic- 
ture of Keplerian orbital motion, a binary system of 
masses should show an inverse correlation between the 
ratio of their velocities, V1/V2, and the ratio of their 
masses, Mi / M^ i .e. V1/V2 = MijM\. As shown in 
IWang et al.1 ((2009) for their sample of double-peaked 
[O III]A5007 AGN, the ratio of the line offsets is equiv- 
alent to the velocity ratio, and the mass ratio can be 
approximated by the luminosity ratio multiplied by a 
factor representing the ratio of the accretion rates, £1,2 ■ 
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This yields the relation Lb/L r = eb, r AV r / AV b between 
the blue and red components. For our sample, we have 
measured the requisite observational properties to per- 
form this same analysis, but with the emission lines 
[Ne V]A3426 and [Ne III]A3869. 

Figure QT] shows AVb/AV r plotted against Lb/L r for 
both [Ne V]A3426 and [Ne III]A3869. The plots include 
both our sample and the low-redshift comparison sample. 
Also over-plotted are the same confirmed/strong dual 
AGN candidates and extended/unresolved NLR AGN 
shown in Figure [10] and discussed in Section 15.31 From 
a Spearman rank correlation test, it is clear that there 
is a strong negative correlation between AVb/AV r and 
Lb/L r , which is in the same sense as expected based 
on the theoretical relation of the binary, Keplerian or- 
bit. Specifically, a strong correlation is seen individually 
in our sample, the comparison sample, and in the com- 
bined sample. For both [Ne V]A3426 and [Ne III]A3869, 
the comparison sample shows slightly stronger correla- 
tions than for o ur sample. 

As argued in iWang et all (|2009D . in a major merger 
one might expect that the SMBHs will be in similar envi- 
ronments and therefore, on average, should have similar 
Eddington ratios (tb,r — 1)- Therefore, we over-plotted 
the theoretical relation for tb,r = 1) which has the fol- 
lowing form in base-10 logarithm space: log 10 [Lb/L r ] — 
log 10 [e 6r ] - \og w [V b /V r ]. For both [Ne V]A3426 and 
[Ne III]A3869 the colored/starred sources are gener- 
ally consistent with the theoretical and best-fit lines, 
and they span a large enough dynamic range to dis- 
play the same general trend as our sample. The best- 
fit linear relations (weighted by the X and Y errors) are 
\og w [L b /L r ] = (0.16±0.02)-(0.75±0.09)xlog 10 [V b /V r ] 
([Ne V]A3426) and log w [L b /L r ] = (0.05 ± 0.02) - (0.51 ± 
0.04) x \og 10 [V b /V r ] ([Ne III]A3869). These relations are 
similar to the theoretical relation for e b ,r — 1? though 
they are technically not consistent when accounting for 
the la errors, with both fits yielding shallower slopes. 
The shallower slopes appear to be caused by some of the 
sources at the low L b / L r end which extend out to larger 
Vb/Vr ratios than the rest of the sample. In this portion 
of the plot, our sources are systematically shifted above 
the theoretical relation. Finally, we have over-plotted 
the theoretical relation expected for a biconical outflow , 
L b /L r = (AV b /AV r ) 3 (developed in IWang et a"Dl2009h . 
This relation is in the nearly the opposite sense as the 
for the binary relation and is far from agreeing with the 
best-fit relations. In Section 17.2.31 we will discuss how 
this correlation and the offset from the theoretical rela- 
tion may be consistent with the presence of both dual 
AGN and outflows in our sample. 

While th e vers ion of this analysis presented in 
IWang et al.1 ((2009) utilized host galaxy redshifts (ob- 
tained through fitting template galaxy spectra), we can 
not obtain that information for our sample since the 
quasar continuum outshines the galaxy starlight. In- 
stead, the velocity offsets used in our analysis here are 
relative to the SDSS redshifts {zsdss) discussed in Sec- 
tion [3] However, we have attempted to investigate the 
dependence of our results on the choice of redshift. Re- 
calling our discussion in Section [31 we also have red- 
shifts based upon Mg IIA2800 which (in the single 
SMBH scenario) should trace the central, active SMBH 
redshift. Therefore, we examined Vb/V r versus Lb/L r 



for the subset of our sample with robust (percent error 
< 0.1%) ZMgii values. We find that this does not intro- 
duce a significant change in the best-fit coefficients for 
[NeV]A3426 or [Ne III]A3869. 

6. RADIO LOUDNESS 

Since 91% (119/131) of our sources are in the FIRST 
footprint (5cr flux limit of ^750 mJy), to determine 
their radio-loudness we have adopted the commonly 
used definition of radio to optical luminosity ratio 1Z — 
L5GHZ/L2500 A, with 1Z = 1 being t he cutoff value for 
the radio- loud classification (|Ho|l2002l) . From our parent 
sample, the fraction of radio-loud SDSS quasars in the 
FIRST footprint is 10%, whereas that fraction is 23% in 
our final double-peaked [Ne V]/[Ne III] sample. For com- 
parison, [Smith^T!!!] (|2010D obtained 9% (parent sample 
at z < 0.8) and 27% (double-peaked [O III]A5007 sam- 
ple), respectively. This preferential selection of radio 
sources over the parent quasar population in both stud- 
ies suggests that the origin of the double-peaks might be 
related to the presence of radio jets in so me sources, as 
in the case of SDSS J151709.20+335324.7 (|Rosario et al.l 
[20H . 

7. INTERPRETATION 

In this section, we synthesize the results of the previ- 
ous sections to highlight the most likely physical scenar- 
ios that produce the line splitting and line offsets in our 
sample of quasars with high-ionization, double-peaked 
narrow emission lines. In particular, we examine the pos- 
sibilities of AGN outflows and dual AGN. 

7.1. Examining the Outflow Hypothesis 

The correlation between line-splitting and line 
width evident in the blue and red systems of both 
[Ne V]A3426 and [Ne III]A3869 (Figure [5]) indicates that 
the mechanisms producing both the line-splittings and 
the line broadening are related, an observation which is 
consistent with the two emission line components orig- 
inating near the same AGN. In general, emission line 
velocity-offsets from the host galaxy or quasar redshift 
are often interpreted as evidence for outflowing photoion- 
ized gas; powerful AGN, including quasars, are known 
to be cap able of driving high -velocity and/or large-scale 
outflows ([Fischer et al.ll20ril ) . We will discuss how prop- 
erties of our sample are consistent with some of the mech- 
anisms known to produce outflows, and how they result 
in stratified NLRs. Furthermore, we will discuss our in- 
terpretation within the context of a proposed outflow and 
stratification geometry shown in Figure 1121 

7.1.1. Mechanisms Producing AGN Outflows 

Two of the commonly proposed mechanisms for driv- 
ing outflows in powerful AGN are radiation pressure from 
the accretion disk and radio jets. Figure [H] shows that the 
velocity-splittings are correlated with the quasar Edding- 
ton ratio, / E dd, for [Ne V]A3426 and for [Ne III]A3869. 
This result is consistent with the notion that outflows 
can be driven by radiation pressure from an accretion 
disk, and that more actively accreting SMBHs will drive 
stronger outflows. This scenario is shown in Figure [T^] 
with the radiation emanating from the acc retion disk 
in the commonly assumed bi-conical shape ([Antonuccil 
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Fig. 12. — Cross-section of the NLR showing our proposed geom- 
etry of the outfiow/ionization stratification scenario and the rela- 
tive origins of [Ne V]A3426 and [Ne III]A3869. The ISM around 
the central AGN is assumed to be spherically symmetric, the ori- 
entation of the outflow axis is 45° from face-on, and we only show 
the NLR out to a radius which includes all of the [Ne V]A3426 and 
[Ne III]A3869 emission. The bi-cone of ionizing radiation (with 
an opening angle of 45°) emanating from the central SMBH is 
indicated by the gray shaded regions. Various individual por- 
tions of the ISM and NLR are indicated by solid circles. Por- 
tions where [Ne V]A3426 is produced but not [Ne III]A3869 are 
shown as solid white circles, portions where both [Ne V]A3426 and 
[Ne III]A3869 are produced are indicated by solid gray circles, and 
portions where [Ne III]A3869 is produced but not [Ne V]A3426 are 
indicated by solid black circles. Portions of the ISM not in the 
ionization cone are indicated by light gray circles. The black ar- 
rows indicate the velocities relative to the central SMBH, and the 
colored waves / arrows represent Dopplcr shifting of the light due to 
the velocity component along the line of sight to the observer. 

Il993f). Our resu lt is a lso consistent with the results 



of IKomossa et aLl (|2008|) who find that their sample of 
NLSls with offset [O III]A5007 lines have relatively large 
Eddington ratios which might be reflective of the ra- 
diation driving the outflow. Finally, that the correla- 
tion with /Edd is stronger for [Ne V]A3426 than for 
[Ne III]A3869 fits with a picture in which the NLR gas 
closest to the central engine is more strongly accelerated 
by the radiation pressure, as indicated in Figure[TOJby the 
larger velocity vectors on material closest to the SMBH. 

However, it is possible that not all of the outflows are 
driven by radiation pressure. In particular, radio jets 
are known to drive outflows in AGN by cntrainmcnt of 
NLR material, and our sample (and the [O III]-selected 
sample) have radio loud fractions of ~25% (Section [5]). 
Therefore, the high fraction of radio loud quasars in our 
sample relative to the parent quasar sample is a strong 
indication that the double-peaked sample we have com- 
piled includes quasars with outflowing components. 

7.1.2. Evidence for Stratified NLRs 

If outflows are a common mechanism for driving line- 
splitting in NLRs, then this should naturally result in a 
stratified NLR since lines of greater LP. and n cr it will be 
preferentially produced nearer to the AGN where they 
are accelerated to higher velocities relative to lines of 
lower LP and n cr it- The relative origins and velocities 
of emission lines in the stratification scenario are illus- 
trated in Figure [121 where [Ne V]A3426 originates closer 



to the SMBH than [Ne III]A3869. We examined this sce- 
nario by looking at the relationship between the velocity- 
offsets of [Ne V]A3426 and [Ne III]A3869 in Section 
15.31 In that analysis, Figure ITOR shows that 60% of the 
sources have a significantly larger blue velocity offset in 
[Ne V]A3426 compared to [Ne III]A3869, suggestive of a 
stratified NLR. We are only seeing the projected velocity- 
offsets, so this percentage might represent a lower limit 
on the number of sources with stratified NLRs when ac- 
counting for random orientations of the outflow axes. We 
note that in Figure [TO] the largest stratifications are seen 
at the largest velocities. This could also be a result of 
projection effects, since orientations which reduce the ra- 
dial velocity components will also reduce the observed 
stratification. However, this trend may also be due, in 
part, to the physical effect of stronger outflows (larger 
velocities) producing larger stratifications. 

For most of the sources in our sample we can only place 
upper limits on the fluxes and velocity offsets of blue 
[O II]A3727 components. This is partly due to the blend- 
ing of the A3726, 3729A doublet (Section |2T2"|). However, 
if the double emission line components are driven by out- 
flows and the NLR is stratified, one might expect that 
[O II] A3727, with a relatively small ionization potential of 
LP. = 13.614 eV and critical density of n crit = 3.4 x 10 3 
cm~ 3 , will be strongest in a portion of the NLR rela- 
tively further from the central AGN. As a result, it will 
not be accelerated to high velocities, resulting in a small 
line-splitting. 

With the outflow axis oriented at some angle interme- 
diate to edge-on or face-on, the redshifted NLR emis- 
sion will consequently be more attenuated than the 
blueshifted emission, as illustrated by the 'obscuring 
ISM' labeled in Figure[TOJ The portions of the NLR with 
the largest line of sight velocity components will be most 
obscured. Conversely, the portions with the smallest line 
of sight velocity components will be the least obscured. 
The result is that, in the presence of such attenuation, 
the most redshifted portion of the emission lines will ob- 
scured, moving the observed position of the red emission 
peak closer to the systemic, i.e. non-Doppler shifted, 
rcdshift. In this case the red component resembles the 
'classical' NLR. This is consistent with observations in 
which offset narrow emission lines in AGN are usually 
blueward of the systemic velocity, indicating that we are 
able to view the outflowing component moving toward 
the observer, while the component moving away from 
the observer is obscured by a larger column of dust. For 
example, from Figure [7] it appears that the mean magni- 
tude of the red component velocity-offsets from the sys- 
temic redshift are generally smaller than those of the blue 
components, consistent with the notion that the red com- 
ponents are less dominated by outflows. 

Note that in Figure[TOJ with sufficient attenuation even 
the NLR emission which is least Doppler-shifted will be 
obscured, resulting an apparent blueshifting of the red 
component, as is occ asionally seen in som e of our sources 
and in other studies ()Spoon & Holt 2009). For example, 
the blue [Ne V]A3426 component may be emitted from 
a portion of the NLR on the observer's side which is 
closest to the central source and moving at the great- 
est velocity (e.g. white dots in Figure [T2l . while the 
red [Ne V]A3426 component is from a portion which is 
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further from the central source (e.g. grey dots in Fig- 
ure [T2l . We did find in Section 15.21 evidence for a mild 
positive correlation between the red line offsets and the 
quasar Eddington ratio. This suggests that, while the 
red system tends to represent the 'classical' NLR, it is 
still effected by the radiation pressure since it must orig- 
inate close enough to the central source where the ioniz- 
ing flux is sufficient. Additionally, Figures HUB and HUD 
show that there is some evidence for stratification of the 
red systems (though much less significant than for the 
blue systems). 

7.2. Implications for Dual AGN at High-Redshift 

It is possible that some of the sources in our sample 
may host two SMBHs following a galaxy merger. In this 
case, the double-peaks may be from two distinct NLRs 
that each accompanies its own active SMBH, or perhaps 
two NLR pe aks are influenced by the orbital motion of 
two SMBHs (jBlecha et al.ll2012[ ). So far there are only a 
handful of known plausible merger remnants hosting two 
AGN at redshifts comparable to our sample: z ~ 0.70 9 
(jGerke et al]l2007h. z ~ 0.78 (IComerford et al.ll2009bh . 
and z ~ 1.175 (jBarrows et"al~ll2012l ). Since galaxy merg- 
ers were more frequent at higher redshifts, we would like 
to investigate the dual AGN scenario for the sources in 
our sample. 

7.2.1. Sources with No Apparent Ionization Stratification 

The two [O III]-selected confirmed dual AGN for 
which we measured double peaks in [Ne V]A3426 and 
[Ne III]A3869 are consistent with no apparent ioniza- 
tion stratification (Figure [TUj). Additionally, there is only 
one candidate dual AGN identified through the double- 
peaked profile of [Ne V]A3426 and [Ne III]A3869, and it 
is als o consistent with no apparent ionization stratifica- 
tion (jBarrows et al.ll2~012D . In these sources, the evidence 
for outflowing NLR material is less compelling, and the 
line-splitting may instead be produced by orbital mo- 
tion of two AGN about each other. Likewise, the 40% 
percent of our sources plotted in Figure [TUJ with no ap- 
parent ionization stratification may include cases where 
the double-peaks are the result of two SMBHs following 
a galaxy merger. Additionally, we see explicit double 
[O II]A3727 peaks for a subset (11) of our sources which 
perhaps suggests that the outflow scenario is less likely 
in these sources since [O II]A3727 should originate at a 
greater distance from the central AGN, as mentioned in 
Section O 

We note, however, that the lack of an apparent 
ionization stratification does not preclude the possibility 
of an outflow, or general gas kinematic origin of the 
double-peaked emission l ines. For exam ple, the two 
extended NLR AGN from lFu etaD (|2012f l have no mea- 
surable ionization stratification but the double emission 
components are known to be produced by the NLR 
around a single AGN based on integral-field spectroscopy 
and high-resolution imaging. Conversely, a stratified 
NLR (or two stratified NLRs) does not preclude the 
presence of two AGN. For example, as discussed in 
Section [5731 our measurement of [Ne V]/[Ne III]A3869 in 
the confirmed dual AGN SDSS J150243.1+111557 shows 
some evidence for a stratification. 



7.2.2. Dual AGN with Large Velocities 

As is evident from Figure [7] and discussed in Section [4] 
the velocity-splittings in our sample are generally larger 
than those from [O III]-selected samples, which might 
tend to select against likely dual AGN candidates at kpc- 
scale separations since they would not be bound to the 
merging galaxy system with such large velocities. Most 
strong dual AGN candidates have AI^s less than 500 
km s" 1 : AV = 1 50 km s" 1 (jComerford et al.l 120094 
ICivano et al.ll2010h . AV = 500 km s" 1 flXu fe Komossal 
HH), AV = 350 km s" 1 ([Comerford et afl[20Tll), an d 
AV = 420 km s" 1 (jMcGurk et alJl20lH iFu et al.ll2012f ). 

However, there are several candid ates with velocities 
> 500 km s" 1 : = 630 km s" 1 (IGerke et al.ll2007h 
and AV = 700 km s" 1 (jBarrows et~atll2012f K Addition- 
ally, the dual AGN hypothesis could be allowed for larger 
AVs if the AGN pa i rs are at small separations. For ex- 
ample, jBlecha[^tlil] (|2012f ) find in their simulations that 
large AVs (> 500 km s _1 ) are often associated with dual 
AGN at sub-kpc separations during pericentric passages. 
For comparison, 23% of our sample have AV < 500 km 
s -1 , 45% have 500 < AV < 800 km s" 1 , and 32% have 
AV > 800 km s _1 , with a maximum of AV^ = 1665 km 
s _1 . Therefore, though a fraction of our sample exhibit 
AT^s higher than expected for dual AGN, 2/3 fall in the 
range expected for either kpc or sub-kpc separation AGN 
pairs. We note that recent numerical simulations suggest 
dual activation of the SMBHs following a galaxy merger 
is m ost likely to occur at separat ions smaller than 1-10 
kpc (jVan W assenh ove et al.ll2012f ). Therefore, under this 
picture of dual activation, the correlation between Ed- 
dington ratio and line-splitting seen in Figure [9] would 
naturally emerge for a sample of dual AGN. 

As discussed in Section [3l for dual AGN with suffi- 
ciently large orbital velocities the broad emission line 
profiles may be significantly broader than expected 
if both components are Type 1 AGN. To test for 
such additional broadening, we have compared the 
Mg IIA2800 FWHMs with those of Fe II, but find no 
evidence for systematically broadened Mg IIA2800 com- 
pared to Fe II. Furthermore, a K-S test does not indicate 
a significant difference between the Mg IIA2800 FWHM 
distribution of our sample and that of the parent sam- 
ple. However, we note that the mean FW H Mm s \\ value 
for our sample (4740 km s _1 ) is slightly larger than that 
of the parent sample (4580 km s _1 ) which is perhaps 
suggestive of additional Mg IIA2800 broadening. 

7.2.3. Testing the Dynamical Argument 



In Section [5.41 we tested a dynamical argument for the 
presence of dual AGN in our sample. The results are 
generally consistent with the theoretical expectation for 
a binary, Keplerian orbit (FigurefTTjl. However, to under- 
stand the extent to which we can interpret this result, we 
must also strongly consider the role of alternative physi- 
cal scenarios in producing such a correlation. 

First, we note that for both [Ne V]A3426 and 
[Ne III]A3869, based on the Spearman rank test, the com- 
parison sample shows slightly stronger correlations than 
for our sample. If this is an indication that the com- 
parison sample shows stronger evidence for dual AGN, it 
would be consistent with the notion that larger velocity- 
splittings are less likely to be associated with dual AGN 
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since our sample has larger velocity-splittings than the 
comparison sample. 

Second, it is worth noting that the coefficients 
for the best-fit linear relations of [Ne V]A3426 and 
[Ne III]A3869 are in disagreement. This is consistent 
with the correlation being produced by outflows (or at 
least some of the sources experiencing outflows) since 
the red component of [Ne III]A3869 would originate at 
a greater distance from the observer (compared to the 
red [Ne V] A3426 component) and therefore be even more 
obscured relative to the blue component (Figure [i~2|) . 
In this case, the Lf,/L r ratio should be even larger for 
[Ne III] A3869, making the slope shallower as observed. 

However, it is possible that if our sample contains some 
combination of outflows and dual AGN then the outflows 
are responsible for deviations from the theoretical binary 
relation. For example, as discussed in Section 15.41 many 
of the sources at the low Lb/L r , high Vb/V r portion of 
Figure [TT1 are offset above the theoretical relation. These 
sources may be more likely to represent outflows since 
they are trending in the same direction as the outflow re- 
lation. Additionally, they have large Vb/V r ratios because 
the red component is near the systemic redshift, consis- 
tent with attenuation of the redshifted outflow compo- 
nent. Lastly, the Lf,/L r ratios are smaller, indicating 
that the red component is stronger, and the blue com- 
ponent is a lower luminosity, extended wing as is often 
seen in outflows and is seen a few of our sources. If these 
deviant sources are most likely to be outflows, then the 
remainder would be more consistent with the theoreti- 
cal relation. Additionally, the remainder would have a 
Li,/L r distri bution consi s tent w ith tb.r = 1, similar to 
the result of I Wang et ail (|2009l) . 

We note that an additional source of scatter in the cor- 
relation could be due to stochastic accretion, such that 
the luminosity ratio does not accurately reflect the true 
mass ratio. This effect could be particularly significant 
when the SMBHs are at larger separations when gas is 
less efficiently funneled to the nuclear regions. Inter- 
estingly, at z = 0.8 — 1.6, the 3" fiber diameter of the 
SDSS spectrograph corresponds to —22-25 kpc, so that 
our sample may contain such early-stage mergers. 

7.2.4. Estimating the Fraction of Dual AGN in Our 
Sample 

While the completeness of our selection process as dis- 
cussed in Section |4] suggests that there is a significant 
number of double-peaked emitters that we have missed, 
especially at small AT^s, we can not correct for the true 
number since we do not know the shape of the underlying 
distribution. However, at the least we can use our [O III]- 
selected comparison sample to estimate what fraction 
of [O III]A5007 double-peaked emitters that we missed 
based on selection through [Ne V ]A3426 or [Ne III1A 3869. 
Of the 57 Type 1 AGN from (iSmith et al.l l2010l) . we 
could reliably measure double [Ne III] A3869 peaks for 
67% (we note that this fraction is consistent with the 
comparis on between our completeness estimates and 
those of iLiu et all (|2010bf ) in Section |4~2]) . Based on 
this fraction, we arrive at a corrected number of 195 
double-peaked [O III]A5007 sources and a double-peaked 
[O III]A5007 AGN fraction of - 0.5% at 0.8 < z < 1.6. 
This fraction is likely to be a lower estimate because 
the double-peaked [Ne V]A3426 and [Ne III]A3869 lines 



of the [O III] -selected sample were not selected in ex- 
actly the same way (i.e. we had a prior knowledge of 
the double-peaked separation). Therefore, the double- 
peaked [O III]A5007 AGN fraction at 0.8 < z < 1.6 
is likely to be larger, potentially making it consist ent 
with the fraction of ~ 1% found bv lLiu et al.l (|2010b| ) at 
z < 0.8. 

To actually estimate the expected fraction of dual AGN 
in our sample we need the true fractions of double-peaked 
AGN and of dual AGN out of all AGN at 0.8 < z < 1.6, 
neither of which are known. However, we may make 
several reasonable assumptions that provide a rough es- 
timate. First, in order to determine the true fraction of 
double-peaked AGN out of all AGN at0.8<z<1.6, we 
need to correct for both our selection incompleteness and 
ra ndom proje c tion ef fects. This fraction was estimated 
by lShen et al.l (|2011al ) in which they determined that the 
fraction of detectable double-peaked [O III]A5007 AGN 
is only 20% of the actual number of AGN with dou- 
ble [O III]A5007 components. Correcting our estimated 
double-peaked [O III]A5007 fraction of 0.5% yields a 
'true' double [O III]A5007 fraction of -2.5% at 0.8 < 
z < 1.6. The influence of inclination and phase angle 
for selection of dual AGN through doubl e-peaked emis- 
sion li ne profiles was also investigated bv lWang fc Zhou! 
(2012) in which they found that, at a phase angle of 
<j> « 50°, we miss at least 50% of all AGN with double 
emission components. We note that this 50% correction 
is a lower limit since it does not account for instrumen- 
tal resolution, and that applying additional corrections 
based on our completeness estimates would likely make 
th e correction estima ted in this manner similar to that 
of IShen et all (|2011aD . Finally, if we take the dua l AGN 
fraction at z = 1.2 to be ~ 0.05% (|Yu et al.ll201ll ). then 
we estimate the fraction of dual AGN out of double- 
peaked AGN at 0.8 < z < 1.6 to be 2% This fraction is 
several times s maller than the res ults of iFu et all ([2012ft 
(4.5-12%) and IShen et all (|20lTal) (~ 10%) which were 
obtained from follow-up observations of double-peaked 
[O III]A5007 AGN. However, the difference can be at- 
tributed, at least in part, to the small expected number 
of dual AGN at z = 1.2 estimated by lYu et all ()2011l ) 
which is due to the redshift evolution of galaxy morphol- 
ogy in their analysis which yields more late-type galaxies 
with smaller initial SMBH masses at higher redshift. 

A direct test of this through follow-up observations 
is therefore crucial in understanding the frequency of 
galaxy mergers at redshifts z > 0.8 and their role in 
AGN triggering. For example, NIR spectroscopy will be 
capable of accessing the redshifted [O III] A5007 emis- 
sion line for sources in our sample, allowing for a direct 
comparison with the z < 0.8 samples of double-peaked 
AGN. This was done with the z = 1. 175 dual AGN 
candidate CXOXBJ142607.6+353351 in iBarrows et all 
((20121) which was initially selected through double- 
peaked [Ne V]/[Ne III] but for which follow-up NIR 
spectroscopy provided access to [O III]A5007. The ad- 
ditional spatial information of [O III]A5007 provided by 
2D longslit spectroscopy would enable one to determine 
if any of these sources are strong dual AGN candidates. 
Follow-up high-resolution imaging, such as radio observa- 
tions, would be capable of resolving the two AGN cores, 
if present. 
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8. CONCLUSIONS 

We have compiled a sample of 131 quasars at z = 
0.8 — 1.6 which show double emission line compo- 
nents in either of the high-ionization narrow lines 
[Ne V]A3426 and [Ne III]A3869. The purpose of this 
search was to identify high-redshift analogs of the double- 
peaked [O III]A5007 sources found in several previous 
studies. Those double-peaked [O III]A5007 sources are 
believed to represent complex gas kinematics, large-scale 
outflows, or in a few cases dual AGN. Given the increased 
frequency of galaxy mergers at higher redshifts and their 
importance in models of galaxy evolution, we have inves- 
tigated these phenomena at higher redshifts using our 
sample, with the following conclusions: 

o There is a clear bias towards selecting double-peaks 
with large velocity-splittings. This bias was made 
apparent by our comparison of the velocity-offsets 
of the blue and red components in our sample to 
those of [O III]-selected samples, and it is corrob- 
orated by the results of our completeness simula- 
tions. This selection bias is not surprising, and 
it is imposed by the relatively weaker intensities 
of [Ne V]A3426 and [Ne III]A3869 compared to 
[O III]A5007. 

o We have found two results suggesting that both 
the blue and red systems are influenced by kine- 
matics in the NLR. First, the line-widths of both 
the blue and red components are strongly corre- 
lated with the line-splittings, suggesting a common 
origin. Second, we find that the individual offsets 
for both the 'blue' and 'red' systems are positively 
correlated with the quasar Eddington ratio, sug- 
gesting that the SMBH accretion rate and therefore 
the radiation pressure is responsible for driving the 
line-offsets in the blueward direction for both line 
components. 

o We find evidence suggesting that the observed kine- 
matics are strongest in the blue systems. This is 
suggested because the blue systems' have larger ve- 
locity shifts from the quasar redshift, those velocity 
offsets show the strongest correlations with the Ed- 
dington ratio, and the blue systems show the high- 
est degree of ionization stratification. This further 
suggests that the red outflowing components are 
generally more obscured. 

o We have found that a significant fraction (~ 23%) 
are radio loud, compared to the 10% radio loud 
fraction of the parent sample. 

Taken together, the previous conclusions paint a pic- 
ture in which the blue systems originate in a portion of 
the NLR much closer to the AGN where they are acceler- 
ated by the accretion disk radiation pressure or radio jets 
to high velocities. This explains the large blueshifts, the 
stronger correlation with Eddington ratio and the pro- 
nounced ionization stratification. The red system origi- 
nates further from the AGN where it is not accelerated 
to the high velocities of the blue system but is never- 
theless close enough so it's bulk velocity offset is also 
influenced by the AGN radiation pressure. This explains 



the smaller velocity offsets, the much weaker correlation 
with Eddington ratio, and the relatively less pronounced 
ionization stratification compared to the blue systems. A 
generalized schematic of this scenario in Figure [T2] The 
enhanced radio loud fraction relative to the parent sam- 
ple also suggests that radio jets may be another mech- 
anism which is capable of accelerating the NLR clouds 
to produce the line offsets. This sample can be used 
to study outflows from luminous AGN at relatively high 
redshifts when AGN feedback may have been an impor- 
tant factor in the growth of massive galaxies. 

There are several interesting results which leave 
open the possibility for dual AGN in the sample. In 
particular, several of our correlations can be thought 
of as consistent with the dual AGN scenario, and even 
suggest that the sample is likely to include dual AGN: 

o The correlation between velocity-splitting and Ed- 
dington ratio, while consistent with the picture 
of radiatively-driven outflows, could plausibly be 
consistent with orbiting SMBH pairs in which en- 
hanced accretion is more likely to occur at smaller 
separations where the SMBH orbital velocities will 
be largest. 

o We have found that a subset of our sample (40%) 
are consistent with no measurable ionization strat- 
ification between [Ne V]A3426 and [Ne III]A3869, 
similar to other dual AGN and strong candidate 
dual AGN. 

o We have found that our sample shows a corre- 
lation between the velocity-offset ratio and the 
luminosity ratio of the blue and red components. 
This correlation is broadly consistent with the 
theoretical expectation for a binary, Keplerian or- 
bit, though our sample seems to be systematically 
offset above the relation. We have shown how this 
deviation could be produced in a sample which 
includes a combination of AGN outflows and dual 
AGN. 

o We have estimated the fraction of dual AGN out of 
double-peaked AGN that we expect at 0.8 < z < 
1.6, finding a fraction (2%) which is smaller than 
that estimated at lower redshifts. However, we cau- 
tion that a significant - and perhaps the primary - 
reason for this lower fraction is the small estimated 
number of high redshift dual AGN that we adopt. 

Follow-up NIR observations to access the 
[O III]A5007 line in our sources would allow for a 
direct comparison with the [O III]A5007 velocity and 
spatial profiles of the z < 0.80 samples, allowing for a 
more robust assessment of the origin of the double-peaks 
in the high-ionization narrow emission lines. Therefore, 
this sample represents an initial step toward extending 
the study of double-peaked emission line AGN to higher 
redshifts. 
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